Macrospin Tunneling and Magnetopolaritons with Nanomechanical Interference 
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We theoretically address the quantum dynamics of a nanomechanical resonator coupled to the 
macrospin of a magnetic nanoparticle by both instanton and perturbative approaches. We demon- 
strate suppression of the tunneling between opposite magnetizations and destruction of magnetopo- 
laritons (coherent magneto-mechanical oscillations) by nanomechanical interference. The predictions 
can be verified experimentally by a molecular magnet attached to a nanomechanical bridge. 
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The first direct observation of quantum behavior of a 
macroscopic mechanical resonator constituting a nano- 
electro-mechanical system (NEMS) has been reported 
recently pQ. This opens a wide range of new possibili- 
ties for testing quantum-mechanical principles on macro- 
scopic objects and has the potential to impact sensor 
technology. NEMS's have also been suggested to operate 
as qubits and memory elements for quantum-information 
processing [2j. Proposals and realizations of two-level 
systems (e.g., superconducting qubits) coupled to me- 
chanical modes [2, 3j allow quantum measurements on 
the mechanical resonator. Here, we study quantum ef- 
fects in a NEMS coupled to a ferromagnetic nanoparticle 
such as a single- molecule magnet (SMM). 

The dynamics of a magnetic order parameter and a 
mechanical resonator are coupled by conservation of an- 
gular momentum |4J. The magnetization dynamics of a 
ferromagnetic particle |4J as well as macrospin tunnel- 
ing oscillations in an SMM [5j should in principle induce 
magnetomechanical motion. However, the semiclassical 
treatment fails when the coupling becomes stronger and 
quantum mechanical effects such as freezing of spin tun- 
neling [14J manifest themselves, as discussed below. 

In this Letter, we consider a torsional nanomechanical 
resonator [see Fig. [IJa)] consisting of a load (e.g., a mag- 
netic nanoparticle or an SMM attached to a paddle) and 
a mechanical link to the base (e.g., a nanotube or a chem- 
ical bond). The projection of the wave function of the 
macrospin on the two lowest energy levels is equivalent 
to a harmonic oscillator coupled to a two- level system [6]. 
The interference effects discussed here can be understood 
by considering a mechanical resonator in the nth excited 
state [see Fig. [TJb)], which has n-\- 1 probability maxima 
at different torsion angles. For example, the first excited 
state can be thought of as a superposition of two wave 
functions peaked at different torsion angles [7J. 

For the first excited state, the instanton path on a 
unit sphere [the middle geodesic in Fig. [TJc)] is split 
into two equivalent ones [geodesies with a scatter in 
Fig. flic)] due to magnetic anisotropics defining the tun- 




Figure 1: (a) A torsional resonator consisting of a beam (e.g., 
nanotube) and a load (e.g., a magnetic nanoparticle or an 
SMM attached to a paddle), (b) Probability of finding the 
resonator at a given angle for the four lowest energy levels, 
(c) The first excited state of the resonator is effectively a su- 
perposition of two states with positive and negative torsion 
angles. The spin-reversal tunnel path on a unit sphere with 
a rigid resonator (middle geodesic) splits into two equivalent 
ones (fuzzy paths represent uncertainty in the tilt of the me- 
chanical resonator), allowing interference. 

neling trajectory (e.g. an easy xy— plane anisotropy). 
The area between the equivalent paths in Fig. [TJc) mul- 
tiplied by the spin S is equal to the difference in the 
phases accumulated by the two paths, which leads to a 
complete suppression of tunneling for a phase difference 
of 7r. Interference effects are therefore observable when 
S^/h/ I x uj r 7r, where I x is the moment of inertia of 
the load and uo r is the natural frequency of the resonator 
(the feasibility of this regime will be discussed below). 
We find that interference is most significant in systems, 
in which there is a certain ratio between the spin and me- 
chanical angular momenta which is analogous to the se- 
lection (parity) rules in large-angle macrospin-tunneling 
[8]. We predict that the effect is rather robust and can 
be observed at experimentally achievable temperatures in 
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state-of-the-art structures. Furthermore, tunneling can 
be suppressed by raising temperature (thus repopulating 
the lowest states but without increasing decoherence) . 
Below, we derive rigorously that quantum-mechanical os- 
cillations of the resonator indeed lead to a suppression of 
macrospin tunneling and destruction of magnetomechan- 
ical modes. Using the instanton approach, we find that 
the coupling of a magnetic particle with an easy-plane 
anisotropy to a mechanical resonator can only lower the 
tunneling rate, thus stabilizing the spin. 

Consider a magnetic nanoparticle that behaves as a 
rigid, spin-5 object and is characterized by the magnetic 
anisotropy energy H A = E(S x ,S y ,S z ), where S x ^ ViZ ) 
are the spin-projection operators (in units of H). Here, 
E(S X , S y , S z ) is the classical magnetic energy correspond- 
ing to an easy x axis and a transverse perturbation (par- 
ity symmetric about the y — z plane) that couples the 
magnet to the torsional motion (e.g. an easy xy— plane 
anisotropy). Our complete Hamiltonian is 



H = H r + e- iS "*H A e iS "* + jhS x B , 



(1) 



where the spin is coupled to a single mechanical mode 
with frequency uo r and Hamiltonian H r = hou r (a)a + 1/2) 
in terms of creation/annihilation operators a) /a. The 
second term in Eq. describes magnetic anisotropy, 
taking into account its orientation with respect to the 
lattice in terms of the torsion angle cp = (ad) + a*d)/2S 
(a = Sy^2h/I x uj r for the system in Fig. [IJ however, 
the analytical solution presented below will also hold for 
an arbitrary complex a). The third term is due to the 
external magnetic field B along the x axis (7 is minus the 
gyromagnetic ratio). We now transform the Hamiltonian 
(Ilk by a unitary transformation e~ lSx ^ to the rest frame: 

Hr = hhj r (a)a - 



1/2) +H A + jHS x B. 



(2) 



Here, a = a — iS x a/2S. We describe tunneling be- 
tween low-lying states of the macrospin by path inte- 
grals in which coherent states are constructed using the 
Heisenberg-Weyl (the resonator) and SU(2) (the spin) 
groups from a state [t,n) by a standard procedure [9]: 
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|t,n> . 



(3) 



Index n here stands for a Fock state with n phonons 
in the mechanical mode; \ refers to a macrospin state 
pointing to the north pole; 9 and are the Euler angles 
defining direction ft = (#, 0) of the macrospin; and the 
complex- valued z = z r + izi parametrizes a generalized 
coherent state of the harmonic oscillator. 

In the large- S limit, the transition amplitude between 
two approximate eigenstates at = 7r/2, = 0, z = 
ia/2 and # = 7r/2, = 7r, z — —ia/2 can be expressed 
through coherent-state path integrals [TO] as 
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where VQ 

h f dt (z r Zi — ZiZ r ) — S(j){l — cos 9) is the 

netic and Se = — J dtE, where E 
huj r [n + 1/2 + (z- is x a/2)(z* + is x a/2)} + E(0, <j>), 
the potential energy contribution to the action in the 
absence of a magnetic field (s x = sin cos (j)). The term 
Sk describes interference effects that can be exposed 
by treating the transition amplitude in Eq. Q by a 
saddle-point approximation. Each saddle-point path 
acquires a phase S A + 2 A' where A is the solid angle 
spanned by the spin paths connected to the north pole 
by geodesies and A' is the area enclosed by the torsional 
trajectories connected to the origin in the complex plane. 
These phases are known to cause interference effects in 
the tunneling of spins [llj. 

We first calculate the tunneling rate in Eq. Q by a 
quasiclassical treatment in imaginary time. The quasi- 
classical equations of motion in real time minimize the 
action in Eq. Q [T2]: 

Zi = —uj r z r , z r = uj r Zi — (auj r /2) cos <j) sin , 
SOsmO = dE/d(j), 50 sin 6 = -dE/dO . (5) 

In order to find the instanton path in the presence of 
the coupling to the mechanical mode, we integrate these 
equations numerically in imaginary time (t = —ir). The 
splitting of the degenerate modes by the tunnel interac- 
tion can be expressed as: 



A n ocCe Re ( 5 o c( )A. 



(6) 



where S§ 1 = iSk + Se is the Wick rotated instanton ac- 
tion of a quasiclassical trajectory found by solving Eqs. 
([5| for 0, 0, Zi, and z r . C describes the interference of 
different spin trajectories under tunneling selection rules 
(e.g. for an easy xy— plane anisotropy C = cos(7r5)) 
[TT] , For an easy-axis/easy-plane anisotropy described by 
E = Ki cos 2 + K 2 sin 2 sin 2 (K x > K 2 > 0), we find 
that all paths are contained between two extremal paths 
denoted in Fig.[2|b) by dotted lines. The tunnel splitting 
changes from exp(— a 2 /2)Ao [rectangular-shaped dotted 
path in Fig. ^b)] to Aq [cosine-shaped dotted path in 
Fig. ^b)] as we go from the limit uo r <C uoq ~ tj 1 (which 
is the focus of this Letter) to the limit uj r ^> ujo ~ tj 1 
where Ao is the tunnel splitting without coupling to the 
resonator, ujo = S^JK\K2 and tj is the uncoupled (mag- 
netic) instanton tunneling time. Such behavior can be 
understood by eliminating and z r from Eqs. ([5| and 
the imaginary-time Lagrangian: 



L = M x (0)0 2 /2 + M 2 i 2 /2 - C7(0, z) . 



(?) 



(4) 



Here Mi(0) = uJq 1 (1/\^X — V^Asin 2 + a cos 0cj r /2cJo), 
M2 = 2/coV, U((j), z) = — a/A^o sm2 — cj r (^ — cos0a/2) 2 , 
A = K2/K1, and w ir/2 has been used. By inspecting 
the motion of a particle with anisotropic mass in the 
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Figure 2: (a) A particle falls in a camel-back shaped potential 
following the quasiclassical path(bold curve), (b) Results of 
a numerical integration of Eqs. (I5l for different frequencies of 
the mechanical resonator. 



potential U (0, z) [Fig.[2ja)], we see that the quasiclassical 
path connecting two potential-energy maxima leads to a 
smaller tunnel splitting because it experiences a higher 
tunneling barrier compared to an uncoupled instanton 
(this is also expected for other magnetic anisotropics). 

We show in the following that the quasiclassical ap- 
proach only works well for n = but does not capture 
nanomechanical interference effects. The latter can be 
obtained by calculating the quantum fluctuations and by 
retaining the second order terms (depending on the Fock 
number n) in the kinetic-energy of the Lagrangian. In 
the remaining part of this Letter, we consider the limit 
tj <C uo~ l in which the kinetic-energy contribution to 
the action suppresses the resonator dynamics. The res- 
onator contribution can then be calculated by taking the 
matrix element between the initial and final states of the 
mechanical subsystem: 
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m) is the Fock states' matrix 
element of the displacement operator and e lSl ^ h corre- 
sponds to the bare instanton contribution to the path 
integral (without the coupling to the resonator). The ex- 
tra factor in the instanton contribution reflects the phases 
accumulated by multiple paths of the macrospin tunnel- 
ing within the laboratory frame in Eq. 0. These paths 
can destructively interfere, suppressing tunneling at spe- 
cific values of a [see Fig.[3^a)], as becomes clear from the 
following analytical expression |13| : 



P -M 2 /2 



(-iar 




L%- m) (\*\ 2 ), (9) 



where n > m (n < m can be obtained by complex conju- 
gation) and L^T™^ is a generalized Laguerre polynomial. 
We can then anticipate that the tunnel splitting should 
be renormalized by the nanoresonator according to: 



A < A , 



(10) 
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Figure 3: (a) Tunnel splittings as a function of the macrospin- 
resonator coupling a for the first three excited states of the 
resonator. The curves show analytical results, while the 
squares are based on the numerical diagonalization of the 
Hamiltonian corresponding to an Fes SMM. (b) Analogous 
plot for tunnel splittings of the magnetopolariton modes cor- 
responding to the Fock states differing by one phonon. (c) 
Lowest energy levels of the Fes SMM coupled to a mechanical 
resonator obtained by numerical diagonalization. Ao/huu r = 
3 x 10" 6 (the energy is offset by -28.239 K). 



where Ao is the tunnel splitting for a bare macrospin. A 
more rigorous derivation also applicable to the formation 
of magnetopolaritons (m ^ n) is discussed below. 

To tackle the resonant coupling between Fock states 
with arbitrary numbers of phonons, we project Hamilto- 
nian ([2| onto the basis formed by the two lowest energy 
states of the Hamiltonian Ha which is possible when the 
transverse perturbations are small (i.e. when Aq <C co>o ~ 
the distance to the third energy state) . We can represent 
the ground and first excited states split by Aq as [14] 



*± = 4= (h/'x) ± iv-x» 



(ii) 



where states |^±x) represent perturbations of |=Lx). The 
leading-order projection procedure, 



h p = l<MIVv>«vl 



(12) 



= ±X 



leads to the Hamiltonian 1 14 : 



H r + <yHSB 
-e i2S ^A /2 



-e- i2S( ^A /2 
H r - jhSB 



(13) 



where we use (^ m | S x \ip n ) = ±S5 mn for states \^± x ). In 
Eq. (13) we treat Aq as a small perturbation. Therefore, 



in the vicinity of the resonant magnetic field correspond- 
ing to the crossing of Fock states m and n, we can further 
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project Eq. (13) onto the states m and n: 



E n 
■CAo/2 



(14) 



arriving at eigenenergies E n 



(E n + £ m )/2 ± 



^/ — Em) 2 + A^ m /2, where the energy of the Fock 
state n (ra) is = hw r (n + 1/2) + jhSB [E rn = 
huj r (m + 1/2) - jhSB] and A mn — |^c mn |Ao. From 
the eigenenergies, we can immediately see that the mag- 
netopolariton and tunnel splittings are given by A mn , 
which is the main result of this Letter. 

Finally, we confirm our analytical results by an exact 
diagonalization of Hamiltonian in the basis that in- 
cludes 30 Fock states for spin S = 10, which corresponds 
to a 630 x 630 matrix. A typical spin Hamiltonian de- 
scribing an SMM reads 

H A = -DS 2 X + E(S 2 Z - S 2 y ) + C(S% + St) , (15) 

where S± = S z ± iS y . The anisotropy constants were 
taken to correspond to an Feg single-molecule magnet 
with D = 0.292 K, E = 0.046 K, and C = -2.9 x 
10 -5 K |8J. Such parameters result in the tunnel split- 
ting Ao/fcs = 4.5 x 10 -8 K in the absence of coupling 
to the resonator. The resonator frequency is chosen as 
hu r /kB = 3 x 10 5 Ao//c# ~ 15 mK where ks is the 
Boltzmann constant. In Fig. |3jc), we plot the calcu- 
lated lowest eigenenergies. The energy (anti-) crossings 
at zero field corresponding to the tunnel splitting are 
shown in Fig. [3](a) as a function of the coupling pa- 
rameter a. We observe a perfect agreement between 
the results of Eq. ( fTpj ) (lines) and the results of the nu- 
merical diagonalization (squares). At specific values of 
a, we observe destructive interference that completely 
quenches the macrospin tunneling. The tunnel splittings 
in Fig.[3^c) at nonzero magnetic fields correspond to mag- 
netopolariton formation and are only possible for a fi- 
nite coupling to the resonator. The magnetopolariton 
splittings reveal interference effects as a function of this 
coupling (see Fig. [3^b) for analytical (lines) as well as 
numerical (squares) results). 

The tunnel splittings can be measured by the Landau- 
Zener method employed in Ref. |8J. The mechanical res- 
onator has to be cooled to temperatures T <C hwr/ks 
(thus higher frequency resonators are preferable); e.g. 
by quantum-optical [15] or conventional cryogenic tech- 
niques (in Ref. [lj hu r /kB ~ 0.1 K). The critical param- 
eter for the observation of magnetopolariton modes and 
interference effects is the macrospin-resonator coupling a. 
We estimate its value for a device that contains an SMM 
strongly absorbed to a paddle of the size 20 x 20 x 10nm 3 , 
with a single- wall carbon nanotube serving as a mechan- 
ical link [see Fig. HTa)]. For a torsional spring constant 
K = 10- 18 N-m [16J, S = 37 (for a Mn 17 SMM [17]) 
and a moment of inertia I x = 10 -36 kg-m 2 , we obtain 
uj r = ^/K/I x ~ 1 GHz and a ~ 0.02. This coupling 



is large enough for the observation of magnetopolariton 
modes but too small for observing interference effects for 
which a has to be comparable to 0.5 [see Fig. |3ja)]. Suf- 
ficiently large values of the coupling can be achieved in a 
Mni2 SMM bridged between leads [18]. For parameters 
similar to the ones used in Ref. [5j, I x ~ 10 -41 kg-m 2 , 
uo r ~ 1 GHz and S = 10, we arrive at a ~ 1.5. The spin- 
resonator coupling can be increased by lowering the tor- 
sional stiffness and the moment of inertia or by increasing 
the spin S. The energy levels in Fig.[3^c) can also be used 
for quantum manipulations and single-phonon control of 
a mechanical resonator. This, however, requires larger 
tunnel splittings in order to overcome decoherence. The 
tunnel splitting can be easily increased by applying the 
magnetic field normal to the anisotropy axis |19| . 

To conclude, we found a quantum-mechanical solution 
for the coupled motion of a macrospin and a mechanical 
resonator. We study tunnel splittings and avoided level 
crossings corresponding to formation of magnetopolari- 
tons, both of which should be detectable by quantum op- 
tical techniques or by studying Landau-Zener transitions 
[8J. In the strong spin-resonator coupling regime, we pre- 
dict suppression of the tunneling of magnetization and 
destruction of magnetopolaritons by interference of the 
spin tunneling paths resulting from the quantum state of 
the resonator. We predict that the magnetism in SMMs 
can be significantly stabilized against quantum fluctu- 
ations by sticking them to mechanical resonators with 
large quantum fluctuations. Results presented here are 
relevant for possible realizations of quantum control of 
magnetization at a single-phonon level. 
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